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Abstract Extremely Compact Stars (ECS) contain trapped null geodesies. When 
such objects enter the evolution period admitting geodetical motion of neutrinos, cer- 
tain part of neutrinos produced in their interior will be trapped influencing their neu- 
trino luminosity and thermal evolution. We study neutrino trapping in the braneworld 
ECS, assuming uniform distribution of neutrino emissivity and massless neutrinos. 
We give the efficiency of the neutrino trapping effects in the framework of the simple 
model of the internal spacetime with uniform distribution of energy density, and ex- 
ternal spacetime described by the Reissner-Nordstrom geometry characterized by the 
braneworld "tidal" parameter b. For b <0 the external spacetime is of the black-hole 
type, while for > the external spacetime can be of both black-hole and naked- 
singularity type. Then the ECS surface radius R can be located also above the unsta- 
ble (outer) photon circular orbit. Such basically new types of the spacetimes strongly 
alter the trapping phenomena as compared to the standard case of b = 0. It is shown 
that the neutrino trapping effects are slightly lowered by the presence of physically 
more plausible case of < 0, as compared to the standard internal Schwarzschild 
spacetime, while they can be magnified by positive tidal charges ifb<l and lowered 
for ^ > 1. However, potential astrophysical relevance of the trapping phenomena is 
strongly enhanced for negative tidal charges enabling a significant enlargement of the 
ECS surface radius to values coherent with recent observations. 

Keywords Neutrino trapping • Braneworlds • Extremely compact stars 

PACS 95.30.Sf • 04.50.Kd 

1 Introduction 

One of the promising approaches to the higher-dimensional gravity theories is repre- 
sented by the braneworld models, where the observable universe is a 3-brane (domain 
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wall) to which the matter fields are confined, while the gravity field enters the extra 
spatial dimensions of size that could strongly exceed the Planck length scale |7|. 
Such models provide an elegant solution to the hierarchy problem of electroweak and 
quantum gravity scales, as these become of the same order TeV) due to large extra 
dimensions and could be well tested on the planned supercollider experiments 1 15|. 
Gravity can be localized near the brane at low energies even with an infinite size 
extra dimension with warped spacetime satisfying the 5D Einstein equations contain- 
ing negative cosmological constant |32| and an arbitrary energy-momentum tensor 
allowed on the brane |37 1. The Randall-Sundrum model gives 4D Einstein gravity in 
low energies and the Newtonian limit appears on the 3-brane with high accuracy. Sig- 
nificant deviations from the Einstein gravity occur at very high energies, in the early 
universe and in vicinity of compact objects as black holes and neutron stars ll24ll . 

Recently, no exact solution of the full 5D Einstein equations is know, but there 
is a variety of astrophysically plausible special solutions of the 4D effective Einstein 
equations constrained to the brane. Such solutions describe black holes with spherical 
symmetry Iffl or with axial symmetry |5 1 and compact objects that could represent 
neutron (quark) stars fT6|. The black hole spacetimes are determined by geometry of 
the spherical Reissner-Nordstrom (R-N) and axial Kerr-Newman (K-N) type where 
the electric charge squared is substituted by the braneworld "tidal charge" param- 
eter representing the tidal (Weyl tensor) effects of the bulk space onto the 4D black 
hole structure. Astrophysically relevant properties of the braneworld black holes were 
studied in a series of papers devoted to both motion of matter in their vicinity Il42l l6l 
[I]|26l|29| and optical phenomena 1 3l[33l[l2l[l 1 1- 

In the simple model of spherically symmetric stars with uniform energy density 
profile a variety of special solutions having asymptotically Schwarzschildian char- 
acter and satisfying the braneworld boundary conditions were found 1 16|. The most 
popular is the one with external spacetime described by the Reissner-Nordstrom ge- 
ometry with the braneworld tidal charge parameter reflecting the tidal effects of the 
bulk and related to the energy density and brane tension. Its properties were exten- 
sively studied both in the weak field limit |10 9| and strong field limit when some 
restrictions on the brane tension were implied from the data of kHz QPOs observed 
in low mass binary systems with a neutron star [2T|. Here we focus our attention to 
the interior of the compact objects, namely to the phenomena related to trapping of 
neutrinos in the so called Extremely Compact Stars (ECS) admitting (by definition) 
existence of trapped null geodesies [44 1. Usually, the braneworld tension is assumed 
positive while the related tidal charge has to be negative 1 16|, but the negative tension 
and related positive tidal charge are not excluded 1 14], so we consider here both posi- 
tive and negative tidal charge of the compact object. Note that exterior of neutron stars 
with positive tidal charges can be described by both black-hole and naked-singularity 
types of the external R-N spacetime (for details see EI]). 

In the standard (b — 0) internal Schwarzschild spacetimes of uniform energy den- 
sity p5 , 39 1 with radius R < 3GM/c^, bound null geodesies exist being concentrated 
around the stable circular null geodesic ll3l [4n . From the behaviour of the effective 
potential of null geodesies in the exterior, vacuum Schwarzschild spacetimes, deter- 
mining the unstable null circular geodesies at the radius rph = 3GM/c^ (see, e.g., 
ESl ). we can conclude that any spherically symmetric, static non-singular interior 
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spacetime with radius R < rph admits existence of bound null geodesies. The trapped 
null geodesies are then governed by both internal and external barriers of the null 
geodesies effective potential. [^The realistic equations of state admitting the existence 
of the extremely compact objects were found and investigated for neutron stars, quark 
stars and Q-stars Il8l l27ll30p4ll20ll3 ll . For spacetimes with a non-zero tidal charge, 
the situation is more complex as we shall demonstrate in the following. 

The existence of bound null geodesies in ECS has interesting astrophysical con- 
sequences. For example, using the notion of the optical reference geometry Il4ll40l 
[38 | it was shown that trapped modes of gravitational waves could influence some 
instabilities in these objects lIJ)- Trapping of neutrinos can be of high importance in 
the interior of the ECS: first, it will suppress the neutrino flow as measured by dis- 
tant observers and, second, it can influence cooling of the ECS in a layer extending 
from some radius depending on details of their structure up to the surface radius. 
The cooling process could even be realized in a "two-temperature" regime, when the 
temperature profile in the interior of the star with no trapped neutrinos differs from 
the profile established in the external layer with trapped neutrinos ||44| modifying 
substantially the standard picture of the neutron star structure as given in II17II18I|T91 
06 45 1. In a different context, trapping of neutrinos inside rotating neutron stars has 
been discussed in ll25l . 

The approximation of free, geodetical motion of neutrinos in the internal space- 
time could be used when the mean free path of neutrinos X > R. Neutrinos have in- 
elastic scatter on electrons (muons) and elastic scatter on neutrons. The scatter cross 
section on electrons (neutrons) ffe (On) determines the mean free path by formula 
X ~ ((7,n,)^' where n,- (/ = e,n) denotes the number density of electrons (neutrons). 
It was shown ll36l that 

while 

A.~3«0'^(i»^)%.. ,2) 

There is > 10 km for Ey < 20 MeV and > 10 km for < 500 keV. There- 
fore, in a few hours old neutron star, see [23 ,361, at temperatures T < 10^ K {Ey ^ 
100 keV), the neutrino motion could be considered geodetical through whole the in- 
ternal spacetime. 

Bound neutrinos with mean free path ^ R will slow down the cooling. Of course, 
they will be re-scattered due to finiteness of the mean free path. An eventual scattering 
of trapped neutrinos will cause change of their impact parameter, therefore, some of 
them will escape the ECS, suppressing thus the slow down of the cooling process in 
the region of neutrino trapping. Clearly, the scattering effect of the trapped neutrinos 
is a complex process deserving sophisticated numerical code based on the Monte 
Carlo method (we expect modelling of this effect in future). 

' In principle, the bound null geodesies could exist also in objects having R > 3GM/c^, e.g., in some 
composite polytropic spheres 1301 . 
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The cooling process deserves sophisticated analytical estimates and detailed nu- 
merical simulations. As a first step in considering the role of trapped neutrinos in 
ECS, efficiency of the neutrino trapping effect was studied in the simple case of the 
internal Schwarzschild spacetime with uniform distribution of energy density (but a 
nontrivial pressure profile) and isotropic and uniform distribution of local neutrino 
luminosity, when all the calculations can be realized in terms of elementary functions 
only (see |i44l).|^The influence on the neutrino luminosity of the star is given by a 
luminosity trapping coefficient relating the total number of trapped neutrinos and the 
total number of radiated neutrinos (per unit time of distant observers). The influence 
on the cooling process is given by two "cooling" trapping coefficients: a "local" one 
given by ratio of trapped and radiated neutrinos at any radius where the trapping oc- 
curs, and the "global" one giving ratio of trapped and radiated neutrinos (per unit time 
of distant observers) integrated over whole the region where the trapping occurs. 

In the present paper we study the role of the braneworld "tidal charge" parameter 
on the existence of ECS and estimate its role in the efficiency of the neutrino trapping 
process using the same trapping coefficients as those introduced in | 44| . We consider 
the trapping effect in the simple case of objects with uniform distribution of energy 
density. The external field is given by the Reissner-Nordstrom geometry fW]. This 
special solution is most extensively studied in the literature, since its external field is 
identical to those of the black hole solution with spherical symmetry lfT4il — there- 
fore, the results of our study could be efficiently compared to the results obtained in 
studies of other phenomena 1 10 21 1. We consider both positive and negative values 
of the braneworld tidal parameter — usually, only the negatively valued tidal charges 
(corresponding to positive tension of the brane) are discussed fT6l|, but the inverse 
situation with positive tidal charge is not excluded and will be treated here, since it 
could demonstrate some new interesting effects. 

Our paper is organized as follows. In Section |2] we summarize properties and 
matching conditions of the internal uniform energy density spacetime and the Reiss- 
ner-Nordstrom external spacetime. In Section |3] null geodesies of both the internal 
and external spacetime are described in terms of properly given effective potential 
and the ECS are classified according to the properties of the trapping region of the 
null geodesic motion. In Section|4] the trapping of neutrinos is discused. In Section|5] 
the trapping efficiency coefficients are defined for both the total neutrino luminosity 
and neutrino cooling process, and determined for all kinds of the braneworld ECS. 
In Section [6] concluding remarks are presented. Throughout the paper, we shall use 
the high-energy units with c = = not stated otherwise. For simplicity, we 
assume zero rest energy of neutrinos, isotropic emission at each radius, and the period 
of evolution of the compact stars, when the temperature is low enough to describe the 
motion of neutrinos by the null geodesies of the spacetime. 



- The internal Schwarzschild spacetime can well represent the spacetime properties of realistic ECS flj] 
I18I19I which are crucial in estimating the role of neutrino trapping. Such a simple spacetime could serve 
quite well as a test bed for realistic models of ECS. 
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2 Braneworld neutron stars 

We describe the braneworld neutron stars by the simple model of uniform energy 
density interior and the "tidally" charged exterior. We give the internal and external 
geometry and the matching conditions. Accordingly, we determine conditions on the 
existence of ECS depending on the tidal parameter characterizing the external space- 
time. 



2. 1 Matching conditions of internal and external spacetimes 

In the standard Schwarzschild coordinates and the high-energy units, the line element 
of the spherically symmetric spacetimes reads 



r^dQ. 



(3) 



The internal solution is characterized by the uniform energy density distribution 
— p — const, and by the tension of the brane — X. We assume p > 0, but we allow 
both possibilities A > 0, A < for completeness. The internal geometry is matched to 
the external geometry at the surface of the star r — R. The line element of the internal 
solution is given by the metric coefficients (r), B (r) that are determined by lfT6l 



A-(r) = 



A{R) 



and 



{B-{r)y 



2GM 



2X) 



(4) 



(5) 



where M = ^npR^ . 



The pressure radial profile is given by 

p{r) ^ [A{r)-A{R)]{\+p/X) 
p 3A{R)-A{r) + [iA{R)-2A{r)]{p/Xy 



(6) 



For r = R, there is p{r)/p = 0. The maximum of the pressure profile is at r = 0. 

The reality condition on the metric coefficient B (r) (taken at r = R) implies a 
relation between A, p and R that can be expressed in the form 



R 



GM ^ p 
-2GM - X' 



(7) 



Considering the restriction R > 2GM (R < 2GM), we can see that the reality 
condition (j?]) is satisfied for all A < (forbidden for all A > 0), while for positive 
tension A > (negative tension A < 0), we obtain a limit on the positive (negative) 
tension given by lfT6l 

fR-2GM\ 

(Notice that considering possibility of R < 2GM, we obtain a limit on the negative 
tension.) 
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The line element of the external geometry is given by the metric coefficients 
A+(r), B^{r) that are determined by 



.2 ' 



(9) 



where, due to the matching conditions on the neutron star surface, i.e., A^{R) ~ 
A^(R), B~{R) = the external mass parameter ^ and external tidal charge 

parameter q are related to the internal geometry parameters p, X (and M) by the 
relations 



-3GMR-, 



:M 1 



(10) 



(11) 



For A > 0, the tidal charge ^ < and < M, while for A < 0, there is <7 > 
and > M. Notice that for A > 0, the condition p < A has to be satisfied in order 
to have ^ > 0. 

For our purposes, it is convenient to express the internal spacetime coefficients us- 
ing the parameters of the external spacetime that can be directly determined from ob- 
servations of accretion and optical phenomena in vicinity of the neutron stars. Since 
the matching conditions imply the relations 



3GMRp 3G.^R 



A 



^.1 
P 



(1-A/p)' 



3G^R 
1 



the internal metric can be expressed in terms of the external parameters 
a new parameter 



X 



3GJ^R 



in the form 



{B-{r)Y=A-\r) = 



1- 



IGJiP- 



R^ 



-1^ 



A-(r) = 



A {R) [34 (R) (2X - 1 ) + 4 (r) (1 - 3X)\ 



2A{R) {2X-l)-A{r)X 
The pressure profile can be expressed in the form 



Pir) 



[A(r)~A{R)]{2X-l) 



(12) 

(13) 
, R, q and 

(14) 

(15) 
(16) 

(17) 



p [3A{R)-A (r)] {X-l) + [3A [R) -2A{r)]X' 

The observational restrictions obtained from the measurements outside the neutron 
star then can be applied to the internal parameters using the relations ([T2|-([T4|. 
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2.2 Limit on existence of uniform density stars 

The limit on the existence of the uniform density spherical configuration is related 
to their compactness and is determined by the condition of pressure finiteness in 
their center We then find the limit on compactness of the star given in terms of the 
gravitational mass related to the internal geometry lfT6l : 



GM 4 

~R -9 



The lowest order correction is given by 

GM 4 

< - 

R - 9 



(l + f)'J 
3p 



1- 



4A 



(18) 



(19) 



Using the external gravitational mass parameter the compactness limit (I81 is 
transformed to the form 



< 



P\2 



:i+f) 



hitroducing the gravitational radius 



rg — G^ 



and dimensionless braneworld tidal charge 



h = 



q 



(20) 



(21) 



(22) 



we can put 



and we arrive to the relations 



X = 



3G^R 3 R 



q 



A^(R) = l-2^+b 
R 



R 



In terms of dimensionless units (rg = 1), the pressure function reads 

P(rRb)- iMr)-A{R)){2b-3R) 
p^''' 3A{R){2b-3R)-3A{r){b-R)' 



(23) 

(24) 
(25) 

(26) 



The pressure increases monotonously with radius decreasing. The central pressure is 
given by 

P , „ „ .^ 3R-2b 



(r = 0,R,b) = 

P 3R(^b+{R-2)R + R-3 



6b 



(27) 



8 



The reality condition of the central pressure reads b+ {R — 2)R > and is equiva- 
lent to the condition A {R) € M. The pressure in the center of the star must be finite 
and positive. These conditions imply the limits ^min(^) on the existence of uniform 
density stars. The central pressure p{r = 0,R,b) diverges when the surface radius R 
satisfies the condition 

3R(^/b+{R-2)R + R-3) +6b = 0. (28) 



This leads to a cubic equation relative to R which gives one solution that is relevant 
for b < 1 and reads 



{b-9){b-l) ^ 




8.y{b-lfb^+b[bib+n)-45] +27 



8W(/7-l)2/73+^[/,(^+17)_45]+27+/7 + 3 . (29) 



The solution is depicted on Figure [T| as the part of ^min(^) denoted by p ^ o°. For 
^7 = 0, we arrive at the standard condition R > (see [44J). For b > I the relevant 
limit is given by the condition of positiveness of the central pressure. It reads 



Rminib\b>l) = h (30) 

and is denoted by p — > 0. For smaller R the central pressure is not positive. 

Note that for 3/4 < ^ < 1 there is a region of surface radii given by the condition 
2b/3 < R < I — y/l—b where p{r = 0) is positive. However, the star surface is lo- 
cated under the inner horizon of the external spacetime that belongs to the black hole 
type R-N spacetime. Such a configuration has to be hidden under the inner horizon 
of the external spacetime and is irrelevant for our considerations. 



2.3 Extremely compact neutron stars 

The ECS (neutron, quark, or hybrid) are defined by the existence of trapped null 
geodesies in the interior of the compact star In the case of internal uniform density 
Schwarzschild spacetimes {b = 0), such objects appear just when the surface of the 
compact star is located under the photon circular geodesic of the external vacuum 
Schwarzschild spacetime located at 3rg |44|. In the braneworld uniform density 
compact stars this simple rule does not hold and the situation is more complex due to 
different character of the external spacetime and its relation to the internal spacetime. 

Assuming the external spacetime to be of the Reissner-Nordstrom (R-N) type 
with the tidal charge substituting the electric charge squared appearing in the standard 
R-N spacetimes, the photon circular geodesies are given by the condition 



(31) 
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However, it is convenient to use dimensionless tidal charge b and dimensionless 
radial coordinate defined by 

r/rg ^ r. (32) 

For negative tidal charges {b < 0) external spacetimes of the black-hole type are 
allowed only ||2T1 . The radii of the photon circular orbits are given by 




However, only the outer photon circular geodesic at rph+ is physically relevant and 
we see immediately that {b < 0) there is rph > 3. 

For braneworld neutron stars with a positive tidal charge (b > 0) both black-hole 
and naked singularity R-N spacetimes are relevant 11211 . In the black-hole type space- 



times (b < 1), the photon circular orbits are given by the relation (33 i, but only the 
outer solution corresponding to an unstable orbit can be astrophysically relevant since 
the surface of the compact star has to be located above radius of the outer horizon. 
In the naked singularity spacetimes (b > 1) two photon circular orbits are can be 
astrophysically relevant when b < 9/8, the inner one being stable, the outer one — 
unstable. In the spacetimes with b >9/di, no photon circular orbits can exist. There- 
fore, in the braneworld spacetimes, the existence of ECS is governed by the interplay 
of the behaviour of the internal and external effective potential of the motion and will 
be determined and classified in the following section. 

3 Null geodesies of internal and external spacetimes 

In terms of the tetrad formalism the metric (|3]l reads 

ds^ = - [a)(')]2 + [a)(''']2 + [©("'l^ + [a)(^)]2, (34) 
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where 

=A(r)df, CO*'"' =B(r)dr, rd0, co**^' = rsin0d(/). (35) 

-1 

, the tetrad of 4-vectors is given by 

^W~AW^' ^^^^ 

Tetrad components of 4-momentum of a test particle or a photon are determined by 
the projections — Py^e^f^^y p^"^ — p^CO^^ which give quantities measured by the 
local observers. 



Since there is e^^j = O)^ 



3. 1 Null geodesies and effective potential 

We consider the period of evolution and cooling of ECS when their temperature falls 
down enough that the motion of neutrinos can be considered free, i.e., geodetical. 
We can assume this period starts at the moment when mean free path of neutrinos 
becomes to be comparable to the radius R, i.e., in hours after the gravitational collapse 
creating the compact object lil7lll8lfT9P3]l36ll . In fact, there are arguments that this 
condition starts to be fulfilled about 50 s after collapse to a proto-neutron star 1 23 , 22] . 
Weak interaction of ultrarelativistic (massless) neutrinos implies their motion along 
null geodesies obeying the equations (A is an affine parameter) 

^=0, p%=0. (37) 

Due to the existence of two Killing vector fields: the temporal d/dt one, and the 
azimuthal d/dcj) one, two conserved components of the 4-momentum must exist: 

E — —pt (energy), L — p^ (axial angular momentum). (38) 

Moreover, the motion plane is central. For a single-particle motion, one can set 9 = 
7t/2 — const, choosing the equatorial plane. 

The motion along null-geodesics is independent of energy (frequency) and can 
conveniently be described in terms of the impact parameter 

i^^. (39) 
E 



Then ( 37 1 yields the relevant equation governing the radial motion in the form 



ipy =A-\r)B-\r)E' ( 1 -A^(r)^ ) . (40) 
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Clearly, the energy E is irrelevant and can be used for rescalling of the affine param- 
eter X. The radial motion is restricted by an effective potential related to the impact 
parameter and defined by the relations 

t<V,n^{ ^2^^ (41) 

^-=(aW 

is the effective potential of the nuU-geodetical motion in the internal braneworld 
spacetime and V^^' is the effective potential of the null-geodetical motion in the ex- 
ternal, vacuum R-N spacetime. We shall see in the following that due to the bulk 
space tidal effects on the matching of the internal and external spacetimes (see [fl6ll 
for details) we obtain a non-standard variety of relations of the effective potentials in 
the ECS interior and exterior 

Using the dimensionless radial coordinate expressed in terms of the gravitational 
radius (rjG./M r), and dimensionless tidal charge h, we obtain the relations 



pff — n 

9\b^{R-l)R\ [{R-b)+R{2b~3R)ZY 



(42) 



^4 



= (43) 

where 

Y = R^A (r) = ^r^{b-2R)+R'^, (44) 

Z = RA (R) = ^/b+{R-2)R. (45) 

Circular null geodesies, located at r^j;) and r(,(e), respectively, are given by the 
local extrema of the effective potential (dVf,ff/dr — 0). In the internal spacetime we 
have to solve a nontrivial equation r. 

d[r^ir,R,b))/dr^ 

[2rR^ {biY - 4-RZ) + 6R^Z) {-b^ {r^{Y - ARZ) + M^Y) + b^Rx 
X (r^ (37 - URZ) + R^ {6R^Z - 9RY + 407) ) + bR^ {R^ {-13R^-Z+ 
+25RY - 667) - 2r^{Y - 6RZ)) + 6R^ {R^Z - 3RY + 67) | x 



9{b+{R- 2)R) Y{R{Y- 3RZ) -biY- 2RZ)f 



0. 



(46) 



This equation can be solved by using numerical methods and determines the loci 
r(.(i) (/?, b) of the internal circular null geodesies that are stable and correspond to a lo- 
cal maximum of the internal effective potential. From the radius r^.t{\ we easily obtain 
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the impact parameter {R,b) using equation (42 1. The existence of the internal sta- 
ble circular null geodesies is allowed only for the surface radius R limited by values 
explicitly given by 

3 



Ri± 



■± 



(47) 



that has to satisfy naturally also the condition R > Rminib)- 

In the external spacetime, the photon circular geodesies and the related local ex- 



trema of the effective potential are given by equation (33 1 and their properties were 
discussed in the previous section. 

When we take into account the (non-)existence of the local extrema of the inter- 
nal/external part of the effective potential Veff, we obtain the whole region of existence 
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of ECS. The function (b) is determined by the relation 



^max-|^^^ ^V_^ /9_|^' for 0</7< 27/20 ' ^^^^ 



rph = 5 (3 + V9^) for b<0 



while for /;ECS(^) 

we have 



^mhf- j^^_^3_^9_|^ forl</,<27/20 ■ ^^'^^ 

The region of ECS in the parameter space is represented in Figure [T| 

Typical behaviour of the effective potential of the nuU-geodetical motion Vgff is 
demonstrated in Figure[2]for appropriatelly chosen values of the tidal charge b. Here, 
and henceforth, we express the radii in units of the gravitational radius rg — G./M . 
The selection of the tidal charge values used in Figure|2]demonstrates the full classi- 
fication of the behavior of the effective potential in both internal and external space- 
times. When the effective potential in the internal or external spacetime (or in both 
of them) has a local extreme corresponding to a photon circular geodesic, trapped 
null geodesies can appear if the surface radius R is properly chosen giving an ECS. 
On the other hand, we find an ordinary compact star if the surface is chosen in such 
a way that no local extrema of the effective potentials exist (see Figure |2^). Sepa- 
ration of the zones of compact stars of different character, both extremely compact 
and ordinary compact, in the parameter space (b — R), is determined by the functions 
^min(^), ^max(^) and ^ph(^) where the last function governs radius of the photon 
circular geodesies in the external spacetime (both of them for 1 < < 9/8). The 
subdivision of the Zone IV in the tidal charge range 1 < < 9/8 is given by the re- 
lation of the magnitude of the effective potential at the surface of the ECS and the 
magnitude of the external effective potential at its local minimum, and is determined 
numerically (see Figure[T]l. 

We divided the region of the parameter space R — b corresponding to the exis- 
tence of ECS into five zones. The classification is based onjhe ejdstence of local 
maxima/minima of V ^'J.J^^'"" in the following way (see Figures 1 andPb 

• Zone I there exist maximum of V^^^ and minimum of V^, both located at r R, 
minVgg < Vgg = R), and there is no local maximum of atr>R; 

• Zone II there exist maximum of V™^^ ai r = R and minimum of ai r ^ R, 
minV^fj < maxy^"}:, and there is no local maximum of V^^\ at r > R; 

• Zone III there exist maximum of at r < R and minimum of V^ff at r — R; 
there is no local maximum of Vl^f at r > R; 

• Zone IV there exist both the local maximum and local minimum of V^f at r > R, 
this zone can be divided into two parts 

a) minimum of y^|' < Vff {r ^ R); 

b) minimum of y^|' > V'^], {r = R). 



R,mn for < 1 
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4 Trapping of neutrinos 

In ECS, some part of produced neutrinos is prevented from escaping these static 
objects because of stable circular null geodesies existing in the internal spacetime, 
and/or because of unstable circular null geodesies existing in the external spacetime 
(Figure|3]l. For external braneworld naked-singularity spacetimes with tidal charge in 
the interval 1 < < 9/8, an additional stable circular null geodesic exists that can be 
important for appearance and structure of accretion discs ([43]), but it is irrelevant 
for trapping of neutrinos radiated by ECS. 

The relation dV^"ydr = determines location of the stable circular null geode- 
sies of the internal spacetime and implies the impact parameter which corresponds to 
the local maximum of the effective potential Vf^ at rj,)-;) to be given by i^^^y 

For < 9/8, the local minimum of V^^ is located at 

'-c(e) = ^(3 + \/9^) (50) 

and corresponds to the unstable circular null geodesies of the Reissner-Nordstrom 
external spacetime, with the impact parameter determined by 



^2 (3 + V9^) 

8(3-2/7 + 79^^) ■ 



''c{e) = W7:. — ^, , /T^^ v- (51) 



For the ECS spacetimes corresponding to all of the Zones I-IVa b of the parameter 
space (b — R), typical behaviour of the effective potentials V^^ and V^^ is represented 
in Figure [3] where the regions of trapped neutrinos are shown explicitly. 



4. 1 Regions of trapping 

Bound neutrinos (depicted by the gray filled areas in Figure |3]l may be fully trapped 
inside the internal spacetime, or may, in some spacetimes, appear outside the interior 
of the ECS being trapped by its strong gravitational field, and entering the interior 
again; note that the trapped neutrinos moving outside the compact star occur in the 
spacetimes of the Zones I, II and IVa. Therefore, we divide the trapped neutrinos into 
two families: 

• "Internal " bound neutrinos (depicted using light gray filled area) with impact pa- 
rameter between ^1[^^{R) and t^^^^y motion of these neutrinos is restricted to the in- 
terior of the ECS. The internal bound neutrinos appear exclusively in spacetimes 
of the Zones III and IVb, but they occur also in the spacetimes of the Zones I and 

IVa. 

• "External " bound neutrinos (depicted by the dark gray filling) with impact pa- 
rameter between i'^^^j and i\JyR)\ such bound neutrinos may leave the ECS inte- 
rior, but they re-enter it. In spacetimes of the Zone II all bound neutrinos are of 
this type. 
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The pericentre and the impact parameter of the "external" marginally bound neu- 
trinos Zones I, II, IVa) can be determined from the condition 

while the pericentre and the impact parameter of the "internal" marginally bound 
neutrinos (ri,(i) ; Zones I, III, IVa, IVb) are given by the condition 

C^('--M0^^^^) = ,;T(^' (53) 

see Figures [3] |4]for the graphical representation. For completeness, we show in Fig.|4] 
also loci r(,(i^ of the stable circular null geodesic and the surface radius R of the ECS. 
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Fig. 4 Location of the radii relevant for the trapping of neutrinos y\y(t), ''b(e) ™d ''c(e) given as functions of 
the surface radius R for characteristic values of the tidal charge b. The range of the radius R corresponding 
to the Zones I-IV of the parameter space {b — R) is depicted in all presented figures. 



4.2 Directional angles 

Considering (without loss of generality, as stated above equation (j39]l) an equatorial 
motion, we can define the directional angle relative to the outward pointed radial 
direction measured in the emitter system (i.e., the local system of static observers in 
the internal spacetime) by the standard relations 

p(^) p{r) 

sinv/=-^, cosv/=-^, (54) 

where 

p(«)=p''4"', P(a)=PM<„) (55) 

are the neutrino momentum component as measured by the static observers. Besides 
conserving components (38 1, and pg — 0, equation (40 1 implies 

/ f-\ 

Pr^±EA-''B{l~A^-^\ . (56) 
For the directional angles we thus obtain relations 



sinv/ = A-, cosVA = ± (l -sin^i/) , (57) 
where A is given by equation (|4]i. 
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Fig. 5 Schematic illustration of the bound-escape ratio at a radius r e {r^^^yR) of an internal spacetime. 
Direction of the neutrino motion with respect to the static observers is related to e^,.) giving the outward 
oriented radial direction. In the spacetimes of the Zones III and IVb only cone is relevant, while for 
spacetime of the Zone II only angle y/e is relevant. 



The directional angle limit for the bound neutrinos is determined for spacetimes 
in the Zones I, II and IVa by the impact parameter ^^^^^ . We arrive to the relation 



cos \j/e = ±- 



1 



^ 9{V9^^ + 3Yz^[{R-b)Y+R{2b-3R)zf ^^^^ 



2\ 2 {-2b + V9^^ + 3) r^R^ [bY + 2R{3R - 2b)zf 

The interval of relevant radii is given by r e {r^i,-^^R). 

The directional angle limit for the "internal" bound neutrinos is determined by 



equation (57i, where for i we use £ = (.mi{R)- In the Zone II spacetimes, there are 



none neutrinos bound only in their interior. On the other hand, for spacetimes in the 
Zones III and IVb, this limit represents the total limit on all bound neutrinos. For the 
internal neutrinos we arrive to the relation 



cos l/i = ± 1 



9/?2 [(R - b)Y+R{2b - 3R)ZY 
P-[bY + 2R{3R-2b)zf 



(59) 



Apparently, the condition y/i > y/e holds at any given radius r < R of the ECS space- 
times belonging to the Zones of the parameter space where both values have good 
meaning. 



4.3 Local escaped to produced neutrinos ratio 



We assume that neutrinos are locally produced by isotropically emitting sources. 
Then escaped-to-produced-neutrinos ratio depends on a geometrical argument only. 
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It is determined by the solid angle 2Q. corresponding to escaping neutrinos (also in- 
ward emitted neutrinos must be involved because even these neutrinos can be radiated 
away), see Figure |5] 

Let A^p, A^e and A^b denote, respectively, the number of produced, escaped and 
trapped neutrinos per unit time of an external static observer at infinity. In order to 
determine the global correction factors 

'^^^'') = A/p(^' ^^^'')-A/p(^ = '-^(^'')' ^^^^ 

it is necessary to introduce the local correction factor for trapped neutrinos at a given 
radius r e ('"b(e);^) (or r € (rbf;),/?), if relevant). 

In the general case of non-isotropic emission of neutrinos, the escaped-to-produ- 
ced ratio is given by the expression 

'ft 2k 

ne{%)= [ //'(f)sinfdfd(/) (61) 







with /7(f) being a directional function of the emission (scattering) process. Be- 
cause of the assumption of isotropic emission of neutrinos in the frame of the static 
observers, the escaping solid angle Qs{%) determines fully the ratio of escaped- 
produced neutrinos and is given by 



fe 2k 



= j j sm^'A^>A^ = 2;r(l -cosfe). (62) 


The escaping correction factor 

, dNJr,R,b) 2Q(wJr.R,b)) 
- d4^ ^ 4. - ^-^osy..W), (63) 

while the complementary factor for trapped neutrinos 

P{r,R,b)^l-e{r,R,b)^ °) ' ' = cos Ye{r,R,B). (64) 

aNp(r,R,b) 

Notice that we consider production and escaping rates at a given radius r, but the 
radius R of the compact object enters the relation as it determines the escaping di- 
rectional angle. The coefficient j5{r,R,b) determines local efficiency of the neutrino 
trapping, i.e., the ratio of the trapped and produced neutrinos at any given radius 
r S (rb(e),^) or r G (ri,(i),/?). Its profile is shown for several representative values 
of the tidal charge b and related sequences of the surface radius R in Figures |6] |7] 
The local maxima of the function l5{r,R,b) (with R and b being fixed) are given by 
the condition dji/dr = which is satisfied at radius r = r(,(i) implying coincidence 
with the radius of the stable circular null geodesic, as anticipated intuitively. For the 
Zone II the maxima are naturally located at the surface radius. In Figures |6] |7] the 
maxima are depicted explicitly. 



19 



"=0. 




0.0 0.2 0.4 0.6 0.8 1.0 
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Fig. 6 Profiles of the local trapping factor P{r/R,b) in the ECS with external spacetime of a black-hole 
type. They are constructed for fixed characteristic values of the brane parameter b and sequences of the 
surface radius R corresponding to the ECS. The upper figure represents the case of negative tidal charges 
(i> = —1.2), whUe the lower figure represents the case of positive tidal charges (b = 0.8). The gray Une 
connects the local maxima of the P{r/R,b) profiles. The P(rjR,b) profiles constructed for parameters 
belonging to different zones of the parameter space are separated by thick lines. In the regions of different 
zones the profiles are constructed for surface radii R spaced by equal distances in the interval of maximal 
and minimal values of R of ECS. 
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Fig. 7 Profiles of the local trapping factor fi(r/R,b) in the ECS with external spacetime of a naked- 
singularity type, constructed for fixed characteristic values of the brane parameter h and sequences of the 
surface radius R corresponding to the ECS. The upper figure corresponds to the case when two circular 
null geodesies of the external spacetime can exist (fe = 1 .08), while the lower figure corresponds only to 
the spacetimes where no circular null geodesic exists in the external spacetime {b = 1.2). The gray line 
connects the local maxima of the J3 (r/R, b) profiles. The figures are constructed in the same way as those 
presented for the black-hole external spacetimes. 
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We can see that the local trapping factor in any given ECS spacetime reaches its 
maximal values at the radius corresponding to the stable circular null geodesic and 
decreases as the radius falls to the centre of the ECS. Nevertheless, the trapping re- 
gion can extend down to the region near the centre only for the most compact ECS 
with external spacetime of the black-hole type (for both negative and positive tidal 
charges), having surface radius close the minimal allowed surface radius. The maxi- 
mum of the j5{r,R,b) profile then reaches value of j3 ^ 1 close to the centre, where 
the stable null circular geodesic is located in such spacetimes. On the other hand, for 
positively tidally charged ECS with external spacetime of the naked-singularity type, 
the j3{r,R,b) profile cannot reach the central part and its maximum is much smaller 
thanj3 = 1. 



4.4 Neutrino production rates 

Generally, the neutrino production is a very complex process depending on detailed 
structure of an extremely compact object. We can express the locally defined neutrino 
production rate in the form 

where d^/K is the number of interactions at radius r, T is the proper time of the static 
observer at the given r, {s/} is the full set of quantities relevant for the production 
rate. We can write that 

d.yr{r) = n{r)r{r)dV{r), (66) 

where n(r), F(r) and dV{r) are the number density of particles entering the neutrino 
production processes, the neutrino production rate and the proper volume element at 
the radius r, respectively. Both n{r) and F(r) are given by detailed structure of the 
extremely compact objects, dV{r) is given by the spacetime geometry. 

Here, considering the uniform energy density braneworld stars (for requirements 
of more realistic model see, e.g., f3jn]45 |), we shall assume the local production rate 
to be proportional to the energy density, i.e., we assume uniform production rate as 
measured by the local static observers; of course, from the point of view of static ob- 
servers at infinity, the production rate will not be distributed uniformly. (According to 
II17II18|[T91 , such toy model could be reasonable good starting point for more realistic 
calculations.) 

In the internal spacetime we can thus write the local neutrino production rate in 
the form 

d^ 

y(r) — p = const. (67) 

dr 

The local neutrino production rate related to the distant static observers is then given 
by the relation including the time-delay factor 

diV 

I = — ^ yA{r,R,b). (68) 
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The number of neutrinos produced at a given radius in a proper volume dV per unit 
time of a distant static observer (that are governed by the internal metric coefficients) 
is given by the relation 

dNp{r,R,b) =I{r,R,b)dV{r) ^4nyA{r,R,b)B{r,R,b)r^dr. (69) 

Integrating through whole the compact object (from to R), we arrive to the global 
neutrino production rate in the form 

f-R 

NJR,b)^AnJ A-{r,R,b)B-{r,R,bydr. (70) 
Jo 

In an analogical way, we can give the expressions for the global rates of escaping and 
trapping of the produced neutrinos: 

f-R 

N^{R.b)^AnJ / {\-co?,\if^{r,R,b))A-{r,R,b)B-{r,R,b)r^dr + 

+A^p('-b(e)), (71) 
fR 

Nb{R,b)^AnJ cos\jfe{r,R,b)A-{r,R,b)B-{r,R,b)r^dr, (72) 



where r^^gjis the radius given by equation (52 1 and cos is determined by 
equation (|58|l. In ECS spacetimes with parameters belonging to the Zones where r^^-^ 
has not good meaning, its role takes rb(i), while we replace cos i/Ze for cos i/i- In such 
ECS spacetimes all bound neutrinos are the internal ones. 



5 Efficiency of neutrino trapping 

In order to characterize the trapping of neutrinos in extremely compact stars, we 
introduce some coefficients giving the efficiency of the trapping effect in connection 
to the total neutrino luminosity and the cooling process in the period of the evolution 
of the star corresponding to the geodetical motion of neutrinos. 



5.1 Trapping coefficient of total neutrino luminosity 

The influence of the trapping effect on the total neutrino luminosity of ECS can be 
appropriately given by the coefficient 3§l relating the number of neutrinos produced 
inside the whole compact star during unit time of distant observers and the number of 
those produced neutrinos that will be captured by the extremely strong gravitational 
field of the star The total luminosity trapping coefficient is therefore given by the 
relation 

J^^^ A-{r,R,b)B- {r,R, b) cos W^{rM dr 



where 



MRM = — 5 , (73) 

^^A-{r,R,b)B-{r,R,bydr 

A-B- ^ ^m{R-b)Y^Ri^lb-m)Z^ 
Y\lR{lb-'iR)Z-bY\ 
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Fig. 8 The total luminosity trapping coefficient ,'^i^{R,h) (full lines) and the internal luminosity trapping 
coefficient JSl (dashed lines) as functions of the surface radius R for fixed characteristic values of 
the tidal charge b. The left plot is given for negative tidal charges — 1.2 < < 0, the right plot for positive 
tidal charges < b < 27/20. The values of b are equally spaced with interval of 0.2 for b < i and 0.7 for 
b>l. 



and the complementary luminosity "escaping" coefficient is determined by the simple 
formula 

<^L{R,b) = l^3§L{R,b). (75) 

We can, moreover, define other global characteristic coefficients. For the "inter- 
nal" neutrinos with motion restricted to the interior of the star, we introduce a coeffi- 
cient 

rR A-/., r, u\r,-^.. D u\ /..E,u\^2, 



l"^^_^A-{r,R,b)B- [r, R, b) cos i/i {r,R,b) dr 



J^A-{r,R,b)B-ir,R,bydr 
and for the "external" neutrinos, we can use a complementary coefficient 



(76) 



(77) 



The results are illustrated for the coefficients ,'^i^{R,b), J2]^{R,b) in Figure |8] 
The behaviour of the total luminosity trapping factor is qualitatively the same for 
ECS with external spacetimes of the black-hole type for both negative and positive 
tidal charges, but it has qualitatively different character for the black-hole and naked- 
singularity ECS spacetimes. In the black-hole-type ECS it monotonously increases 
with decreasing surface radius for any fixed tidal charge (both negative and positive), 
while in the naked-singularity-type ECS spacetimes it reaches a maximum between 
the centre and the surface radius R, and then it decreases to zero value for surface 
radius approaching the minimal value of R. 

For the black-hole-type ECS with negative tidal charges the maximum of the co- 
efficient ^]^{R,b) slightly decreases with decreasing b, while for positively tidally 
charged ECS it increases with increasing b approaching ■^]^{R,b) = I for b-^ I. We 
observe a similar behavior also for J2]^{R,b), but the magnitude of this coefficient is 
smaller in comparison with ^^i^{R,b). In the naked-singularity-type ECS (with posi- 
tive tidal charges) the local maximum of ^]^{R,b) (and ^]^{R,b)) strongly decreases 
with b increasing; of course, for ECS in the Zone III of the space of parameters, the 
parameter ^i^{R,b) is not defined. 
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Fig. 9 The global cooling coefficient related to trapping of neutrinos in the active zone of trapping 
SSQ{R,b). For tidal charge b fixed for values with the same interval as in the previous figure and the 
full range of R corresponding to the ECS is scanned. 



5.2 Trapping coefficient of neutrino cooling process 



The efficiency of the influence of neutrino trapping on the cooling process is most 
effectively described by the local coefficient of trapping relating the trapped and 
produced neutrinos at a given radius of the star that is defined by the relation 



b,{rM = ^{r,R,b). 



(78) 



The local cooling coefficient is therefore given in Figures [6] |7]for appropriately cho- 
sen tidal charges b and related sequences of the surface radius R. All the properties of 
the local trapping coefficient discussed above are thus relevant for the local cooling 
phenomena. 

Further, the cooling process can be appropriately described in a complementary 
manner by a global coefficient for trapping, restricted to the "active" zone, where the 
trapping of neutrinos occurs. The global cooling coefficient is thus defined by the 
relation 



^M,{R,b) 



.A-{r,R,b)B- (r, R, b) cos v/e,i (r, R,b)dr 



,^;i^^^A-{r,R,b)B-{r,R,byAr 



(79) 



The global "cooling" coefficient £i§t^{R,b) is presented in Figure|9] — for few charac- 
teristic values of the tidal charge b, both negative and positive, the full range of the 
ESC surface radius is scanned. 

The behaviour of the global cooling factor .'^^{R,b) is of the same character as 
for the total luminosity factor ■S§i^{R.,b) in both black-hole and naked-singularity-type 
ECS and for both negative and positive tidal charges. Its magnitude is nearly twice 
higher than those of the total luminosity factor for the surface radii R close to their 
maximal value, but it is nearly equal to the total luminosity factor as R approaches its 
minimal value. 
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6 Conclusions 

We have demonstrated that the braneworld ECS can be separated into five classes 
corresponding to five Zones of the parameter space {b — R) where trapping of neu- 
trinos and their motion in the external field of the ECS are realized in qualitatively 
different way. There are fundamental differences of the trapping effect in the field 
of ECS with negative and positive tidal charges. The differences are related to both 
the surface radius of the ECS and the trapping coeffients of both kind — for the to- 
tal luminosity and the cooling process. Moreover, there is a strong difference of the 
trapping efficiency in the ECS with positive tidal charge when related to the black- 
hole spacetimes and the naked-singularity spacetimes. In the first case, the trapping is 
very efficient and the global trapping coefficient can approach the value of 1 for the 
spacetimes with — > 1, while in the second case the coefficients strongly decrease 
for tidal charges growing above the critical value of b — 1 . 

There is an astrophysically very important effect related to the surface radius of 
the ECS with the negative and positive tidal charges. For b <0, the surface radius of 
ECS can significantly overcome the limiting value of R = 3 valid for b ~0. On the 
other hand, for b > 0, the surface radius of ECS decreases under the standard value 
of /? = 3 — for ECS with b ^ 1, there is Rmm ^ 2.3 that is substantially lower than 
the values indicated by observations, and for b > I the ECS surface radius decreases 
further 

In the ECS spacetimes with negative charge there is a crucial effect related to 
the increase of the surface radius of the ECS allowed for decreasing b <0. For tidal 
charges slightly under the value of ~ — 1, the ECS surface radius approaches R — 4 
strongly overcoming the limit of /? 3 corresponding to the standard Schwarzschild 
spacetimes and enters the region of values commonly accepted for observed neu- 
tron stars that can be lowered down to 7? = 3.5 12311131 . In such cases the coeffi- 
cient of total luminosity can be of quite significant value for astrophysically real- 
istic radii of ECS. For example, in the case of ECS with R — 3.5 and b = —1.2, 
there is ^]^{R,b) ^ 0.023. The global cooling coefficient can approach even the 
value of 3§c{R,b) ^ Q.Qll under such conditions, indicating thus very strong and 
significant effects of trapping for braneworld negatively charged ECS, significantly 
overcoming the related values corresponding to the standard (tidally non-charged) 
spacetimes i44i . Because the effect of trapping of neutrinos is cumulative, one can 
expect its relevance in realistic models of ECS to be strongly enhanced by negative 
tidal charges. Note that the conditions on the tidal charge magnitude implied by the 
high-frequency quasiperiodic oscillations observed in some low-mass-X-ray- bina- 
ries containing neutron or quark stars put the limit b < 1.2 (assuming their mass 
m ^ 2Mq) fJTl. This limit allows for sufficient increasing of the ECS surface radius 
with negative tidal charges enabling strong effects of neutrino trapping in the inte- 
rior of such ECS with observationally acceptable radius. The expected lowering of 
the neutrino luminosity in observed neutron stars can be about ^10% and could be 
observationally tested. 

We can conclude that the negative tidal charges cause slight decrease of the ef- 
ficiency of the trapping phenomena in the field of ECS, but they significantly shift 
these phenomena to the region of astrophysically relevant situations due to the strong 
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increase of the surface radius of such ECS, shifting the radius to values acceptable 
observationally. On the other hand, the positive tidal charges generally enhance the 
trapping effects for the ECS spacetimes of the black-hole type, but the surface radius 
of such ECS is strongly shifted to the regions excluded observationally. 
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